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ABSTRACT 



The explosion mechanism behind Type la supernovae is a matter of continuing debate. The diverse attempts to 
identify or at least constrain the physical processes involved in the explosion have been only partially successful 
so far. In this paper we propose to use the thermal X-ray emission from young supernova remnants originated in 
Type la events to extract relevant information concerning the explosions themselves. We have produced a grid 
of thermonuclear supernova models representative of the paradigms currently under debate: pure deflagrations, 
delayed detonations, pulsating delayed detonations and sub-Chandrasekhar explosions, using their density and 
chemical composition profiles to simulate the interaction with the surrounding ambient medium and the ensuing 
plasma heating, non-equilibrium ionization and thermal X-ray emission of the ejecta. Key observational parame- 
ters such as electron temperatures, emission measures and ionization time scales are presented and discussed. We 
find that not only is it possible to identify the explosion mechanism from the spectra of young Type la Supernova 
Remnants, it is in fact necessary to take the detailed ejecta structure into account if such spectra are to be modeled 
in a self-consistent way. Neither element line flux ratios nor element emission measures are good estimates of the 
true ratios of ejected masses, with differences of as much as two or three orders of magnitude for a given model. 
Comparison with observations of the Tycho SNR suggests a delayed detonation as the most probable explosion 
mechanism. Line strengths, line ratios, and the centroid of the Fe Ka line are reasonably well reproduced by a 
model of this kind. 

Subject headings: hydrodynamics — nuclear reactions, nucleosynthesis, abundances — supernovae: general — 
ISM: individual (SN 1572) — supernova remnants — X-rays: ISM 



1. INTRODUCTION 

Type la supernovae (SNIa) play an important role in the 
chemical evolution of galaxies, originating most of the Fe group 
elements in the interstellar medium, and they have become a 
key tool in our understanding of the universe by providing ev- 
idence for its accelerated expansion (Perlmutter et al. 1999). 
Still, our knowledge of the physical processes involved in the 
actual explosions is far from being complete. Fundamental is- 
sues such as the mass of the white dwarf (WD) at the moment of 
the explosion, the location of the ignition and the propagation 
mode of the burning front have not been established yet, and a 
number of different models or paradigms are capable of repro- 
ducing with reasonable accuracy the optical spectra and light 
curves of observed SNIa and their fundamental physical prop- 
erties (Branch and Khokhlov 1995; Hillebrandt and Niemeyer 
2000, for reviews). 

Given this situation, it is important to explore all the po- 
tential sources of information that can shed some light on the 
detailed workings of thermonuclear supernovae. The thermal 
X-ray emission from the shocked ejecta in young supernova 
remnants (SNRs) originated by these explosions is, in princi- 
ple, well suited for this purpose: as the reverse shock advances 
into the ejecta, it compresses and heats them to X-ray emit- 
ting temperatures, producing spectra that depend on the com- 
position, density, temperature and ionization state of the mate- 
rial. The calculation of the expected thermal X-ray emission 
from the shocked ejecta synthesized in the different explosion 
paradigms, however, is not straightforward. The propagation of 



the shocks and the density of the shocked gas is governed by the 
dynamic evolution of the SNR, which at early stages depends 
strongly on the ejecta density profile and is not well represented 
by similarity or unified solutions (Truelove and McKee 1999). 
The densities involved are so low that the ensuing shocks are 
collisionless, and their physics, and particularly how the post- 
shock internal energy is distributed among ions and electrons, is 
not well understood. Moreover, the ages of young SNRs (up to 
a few thousand years) are short compared with the time scales 
for thermal equilibration between ions and electrons and for 
the onset of collisional ionization equilibrium, and therefore the 
plasma is in a transient ionizing state that is difficult to calcu- 
late. 

Our goal is to provide observers with the means to make 
meaningful analysis of the thermal X-ray spectra of young, 
ejecta-dominated type la SNRs and eventually to constrain the 
nature of the event that originated them. Thus, here we present 
the first detailed models for these remnants, evolved from a set 
of theoretical calculations of thermonuclear supernova explo- 
sions. We have focused our attention on the general properties 
of each explosion model, and their consequences on the emit- 
ted X-ray spectra. A detailed analysis of particular objects is a 
delicate task that requires fitting many observables and an ex- 
haustive exploration of the model parameter space in each case. 
This task is left for forthcoming publications. In section 2 we 
give the relevant characteristics of the explosion models in our 
grid, which includes all the paradigms currently under discus- 
sion. The SNR simulation scheme is described in section 3: 
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hydrodynamics, plasma ionization and heating, and synthetic 
spectra. Section 4 contains a discussion of our results and the 
relevant observable quantities that can be derived from the mod- 
els, as well as a preliminary comparison with observations for 
Tycho SNR. The conclusions of this study are presented in sec- 
tion 5. 

2. SUPERNOVA EXPLOSION MODELS 

In order to be able to discriminate between different explo- 
sion mechanisms it is important that all the models be calcu- 
lated consistently, i.e. with the same physics included in the 
same hydro and nucleosynthetic codes. With this requirement 
in mind, we have computed a grid of thermonuclear supernova 
explosion models which is representative of the whole diversity 
of theoretical models currently under debate. All the calcula- 
tions have been performed in one dimension, assuming spheri- 
cal symmetry. For each model, we have followed the explosion 
with a supernova hydrocode until 10^ s after the ignition, when 
the expansion has reached a nearly homologous (r cx v every- 
where) stage. The deposition of the energy of radioactive de- 
cay of ^^Ni on the ejecta is taken into account in this supernova 
code. After 10^ s, however, the ^^Ni disintegration becomes dy- 
namically irrelevant because: 1) most of it has already decayed 
to ^*Co, and 2) an increasing fraction of the energy of the pho- 
tons escapes the supernova due to the drop in opacity caused by 
the expansion. In addition to the detailed nucleosynthesis, we 
have computed the light curves in order to be able to compare 
with historical supernovae. The codes used in the explosion 
calculations are the same as in Bravo et al. (1996). More details 
are given in the appendix. 

The different categories of SNIa explosion models that we 
have included in our grid are: 

• DET: Pure detonation model. In this class of models, 
the flame starts close to the center of the WD, and prop- 
agates supersonically nearly through the whole white 
dwarf, incinerating most of it to Fe-group elements. 

• SCH: Sub-Chandrasekhar mass model. A He detona- 
tion is started at the edge of a He envelope, which feeds 
a converging shock wave into the C-O core. Close to 
the center, the converging shock wave transforms into 
a C-O detonation which propagates outwards and pro- 
cesses the rest of the core. A sandwiched structure is 
produced, rich in Fe-group elements both in the inner 
and in the outer parts of the ejecta (in our model, below 
a Lagrangian mass of ~ O.4M0 and above ~ O.SMq, 
in what was the He envelope, see Fig. Al), and rich in 
intermediate-mass elements plus C-O in between. 

• DBF: Pure deflagration models. In these models the de- 
flagration propagates at the laminar flame velocity (a 
small fraction of the sound velocity) close to the center, 
until the Rayleigh-Taylor instability develops, deforms 
the flame surface, and accelerates the combustion. The 
flame velocity remains subsonic all the way, but when 
the material has reached velocities that are of the same 
order as that of the flame the expansion quenches the 
flame. A large mass of unbumt C-O is ejected in the 
outer layers.^ 

' Other deflagration models can be found in the literature in which the mass of unbumt C-O is small while the mass of intermediate-mass elements is large. This is 
the case, in particular, of the popular W7 model (Nomoto, Thielemann, and Yokoi 1984). 



• DDT: Delayed detonation models. In these models the 
flame propagates initially as a slow deflagration, but 
a transition to a detonation is induced at a prescribed 
flame density during the expansion of the white dwarf. 
As a result, the otherwise unburnt C-O is processed 
partly into Fe-group and partly into intermediate-mass 
elements. 

• PDD: Pulsating delayed detonation models. They differ 
conceptually from delayed detonation models in that the 
transition to detonation is induced only after the white 
dwarf has pulsed. The pulsation is due to the inefficient 
burning produced by a slow deflagration, which is un- 
able to rise the whole energy of the white dwarf (gravi- 
tational + internal + kinetic) above zero. 

Our model grid incorporates examples of aU of these 
paradigms, and it aims to cover the whole parameter space. For 
this paper we have selected eight models from this grid as a rep- 
resentative sample, including extreme cases of DBF, DDT, and 
PDD. Their properties are given in table Al. There, £k is the 
kinetic energy of the ejecta, Mmax and AAI15 are, respectively, 
the bolometric magnitude of the supernova at Ught curve max- 
imum and the change in bolometric magnitude between maxi- 
mum and 15 days later, and the rest of the columns are clearly 
defined by the header. The meaning of the parameters in the 
second column is given in the appendix. 

In Figure Al we show the chemical structure of each model 
after the short lifetime radioactive isotopes have completed 
their decays (all the Fe that appears on the plots, for instance, 
was synthesized as ^^Ni in the explosion). Their density struc- 
tures are shown in Figure A2. In the delayed detonation models, 
the transition from deflagration to detonation happened at a La- 
grangian mass of ~ O.2M0, where its imprint on the density 
profile can be seen. The chemical structure is dominated by 
Fe-group elements up to ^ I.OMq for DDTa, at which point 
the flame density was too low to incinerate matter to nuclear 
statistical equilibrium (NSB), leaving a buffer of intermediate- 
mass elements. For DDTe this buffer is much larger, because of 
the lower densities achieved by the detonation. The outermost 
zone dominated by C-O is also larger than in model DDTa. The 
chemical structure of the pulsating delayed detonation model 
PDDa is very similar to that of DDTa, with the transition to det- 
onation happening at ~ 0.3 M©. The different hydrodynamical 
histories of both models, however, are reflected in their density 
profiles (see Fig. A2). At the time of the transition, the density 
of the external layers of PDDa was on average about two orders 
of magnitude lower than in DDTa. As a result, PDDa displays 
a tail of very low density but larger radii than DDTa. In model 
PDDe the transition to detonation took place at ~ 0.3 M©. Its 
chemical profile is similar to that of DDTe, while its density 
profile is similar to that of PDDa. In the deflagration model 
DEFa, the flame quenched due to the expansion of the WD at a 
lagrangian mass of ~ OJMq, after which a narrow (~ 0. 1 Mq 
wide) region rich in intermediate-mass elements was formed. 
The density structure is very different from the delayed deto- 
nation models, due to the sudden termination of nuclear energy 
generation which results in the formation of a bump of unburnt 
material just above the quenching flame front. In model DEFf 
the results were similar to DEFa, but the flame was quenched at 
~O.9M0. 
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3. SIMULATION SCHEME 

3.1. Overview 

The density profiles from the explosion models described in 
the previous section are used as an input for a standard hydro- 
dynamics code (see section 3.2 for details) that follows the in- 
teraction of the ejecta with a uniform ambient medium (AM). 
Hydrodynamical simulations give the evolution of the radius, 
velocity, density and internal energy per unit mass for each 
fluid element as a function of time, i.e., its dynamic history. 
This dynamic history, together with the chemical composition 
of the fluid element obtained from the explosion models is then 
used to calculate the evolution of the plasma heating and the 
nonequilibrium ionization (NEI) in a self-consistent way (sec- 
tion 3.3), resulting in electron temperatures and ion fractions for 
all the ions of the relevant chemical elements at each time. The 
density, electron temperature and ionization structure within 
a certain region of the SNR at a certain time then lead to a 
synthetic or predicted thermal X-ray spectrum, which can be 
readily calculated with a spectral code (section 3.4) and con- 
volved with any particular instrumental response in order to 
compare the models with observations. In section 3.5 we will 
review the approximations made at each stage of our simula- 
tion method and ascertain its validity for the study of young, 
ejecta-dominated Type la SNRs. 

3.2. Hydrodynamics 

We have built a ID Lagrangian hydrodynamics code such 

as described by Truelove and McKee (1999), with an ideal 
gas equation of state (7 = 5 /3), nonlinear pseudoviscosity and 
no external energy sources or sinks, so the SNR model re- 
mains adiabatic throughout its evolution. We have followed 
the interaction of the ejecta from various explosion models 
with a constant density ambient medium (AM), with p^m = 
10"^'*, 5 X 10"^^ and 10"^^ gcm"^ as a sample representative of 
the interstellar medium (ISM) conditions in most of the Galaxy 
and the Magellanic Clouds. 

3.3. Plasma model and ionization calculations 

The theory for modeling transient heavy element plasmas in 
young SNRs was developed by Hamilton and Sarazin (1984, 
henceforth HS84). These authors describe the plasma with 
three Maxwellian populations: one for the ions (including neu- 
trals), and two for the electrons, a 'cold' component for the 
electrons produced in the ionization process and a 'hot' compo- 
nent for the electrons already present in the preshock gas which 
were heated at the coUisionless shock front. The relative impor- 
tance of the hot electron population thus depends on the amount 
of coUisionless heating at the shock, on the preshock ionization 
state of the ejecta and on the composition of each layer, and 
in general it will decrease with time as more 'cold' electrons 
are produced by postshock ionization, especially if coUisionless 
electron heating is not efficient. The efficiency of this heating 
remains a controversial issue: while theoretical calculations of 
the effect of plasma instabilities on the electron population sug- 
gest a high degree of coUisionless heating at the shock (Cargill 
and Papadopoulos 1988), this is difficult to reconcile with ob- 
servations of high velocity shocks, and this prediction is being 
revised (see Laming 2000, for a discussion of coUisionless elec- 
tron heating). We aUow only for a small amount of coUisionless 
heating at the reverse shock, so we have represented the elec- 
trons by a single MaxweUian population T^. In this two-fluid 



model, the postshock temperature ratio between ions and elec- 
trons, Te^s/Ti,s, is a measure of the efficiency of coUisionless 
heating, varying between (no heating) and 1 (total tempera- 
ture equilibration). For the models presented in this paper, we 
have assumed T^.s/T{,s = 0, but we will explore the effect of dif- 
ferent assumptions in future work. 

All the elements in the unshocked ejecta are assumed to be in 
the singly ionized state. A low preshock ionization state agrees 
with the observational evidence provided by Wu et al. (1983) 
for SN 1006, one of the prototype galactic Type la SNRs, but 
we have found that varying the charge of the unshocked ions 
has Uttle impact in the postshock ionization history, at least 
for charge states below ten. For a more in-depth discussion 
of preshock ionization, see HS84. 

In the shocked plasma, ions and electrons interact 
through Coulomb coUisions which gradually equilibrate the 
Maxwellian populations. At the same time, the ongoing ioniza- 
tion of the ejecta increases the electron population and therefore 
modifies the collision rate. Let p and e be the time-dependent 
mass density and specific internal energy per unit mass of a 
given fluid element as calculated with the hydrodynamic code, 
and He, "i the number densities of electrons and ions, respec- 
tively. At any given time, the total internal energy is distributed 
among ions and electrons so that Ci + ee = £, and the respective 
temperatures are related to the internal energies per unit mass 
in each population by £e,i = OkTe,ine,i)/2p. We can represent 
the chemical composition and ionization state of the fluid ele- 
ment with the normalized number abundances for each element 
X = nx/ni and each ion q fxq = nxi/nx with q = Ofov neu- 
tral and q=Zx for bare ions, so that ^xfx = 1> Y^qfxi = 1- 
The plasma then evolves according to the following equations 
(adapted from HS84): 

_deile _ de^/e _ 1 2^I^TT^I^e\^n-^^\nk in-T^) 
dt dt pe m^mji^l^ / ji_ 

(1) 

= =r~' +Rx'i+^fx'i+^~iIxi+Rxi)fx'i\ , (2) 

dt Amu 

where A = 'Yl,xf^^^ average mass number, Z = 

Ylix ifx'' ^® average ion charge in the fluid element, m^ 
the atomic mass unit. In A the Coulomb logarithm, all funda- 
mental constants take their usual values and all time derivatives 
are Lagrangian. Ix'i and Rxi represent the ionization and re- 
combination rates from ion q of element X, respectively, and 
are taken from the recent compilation by Mazzotta et al. (1998). 
Both equations are obviously coupled in transient plasmas dom- 
inated by heavy elements (where dZ/dt ^G), and they have to 
be integrated simultaneously. At the same time, the processes 
are controlled by the precalculated variation of the hydrody- 
namic quantities p and e. The resulting stiff differential equa- 
tion system is solved with an implicit scheme, resulting in a 
time series of temperatures and ion fractions for each fluid ele- 
ment. The values of Tf. and Ji are recomputed from the updated 
values of Ee,; and ne,i after each new iteration of equations (1) 
and (2). All the chemical elements with fx > 10"^ have been 
included in the calculations. 

3.4. Synthetic spectra 

To produce the synthetic spectra, we have used an updated 
version of the HamUton & Sarazin (HS) code incorporated in 
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the XSPEC software package (Arnaud 1996), which is suitable 
for modeling plasmas that depart significantly from the coUi- 
sional ionization equiUbrium (CIE). Given the density, ioniza- 
tion state and temperature of a fluid element at a certain time, 
the HS code produces a predicted spectrum that can be con- 
volved with the response matrix of any appropriate instrument. 
It is important to mention that this code relies on modem atomic 
calculations for the Fe L shell emission from Li-Uke to Ne-like 
Fe ions (Liedahl, Osterheld, and Goldestein 1995). This code 
includes no atomic data for Ar, and has other limitations which 
restrict its apphcability, but problems of this kind are cormnon 
to most of the currently available NEI codes. The Umitations 
of the spectral code should be considered when comparing the 
predicted spectra with observations, but the quality of our cal- 
culated spectra is satisfactory for our present purposes. 

We will only discuss the integrated ejecta spectra: detailed 
models for the shocked AM spectra and spatially resolved spec- 
troscopy of the ejecta are beyond the scope of this paper, even 
though they can certainly be addressed using the simulation 
scheme we have just described. We have used the response 
matrix for the CCD MOS-1 EPIC camera onboard the XMM- 
Newton observatory, which provides an excellent spectral reso- 
lution and is a reasonable choice for observing extended objects 
like SNRs. The Doppler broadening of fines due to thermal and 
bulk motions have not been included because their effects are 
small at the spectral resolution of XMM EPIC-MOS. 

3.5. Simulation scheme discussion- 
It is important to emphasize that our simulations assume 
spherical symmetry for both the SNIa ejecta and the interaction 
with the AM. Our aim is to show that young SNRs can be a 
valuable tool for studying SNla explosions, and the assumption 
of one-dimensional dynamics is a first step towards that goal. In 
fact, it can be a satisfactory approach in many individual cases, 
especiaUy since Type la SNRs tend to be, on average, more 
symmetric than those originated by core collapse supernovae. 
In contrast to core collapse SNRs, they seem to retain the origi- 
nal ejecta stratification (compare the results of Decourchelle et. 
al. 2001 on Tycho and Lewis et. al. 2003 on N103B with those 
of Hughes et. al. 2000 on Cas A). Whenever deviations from 
spherical symmetry affect the shape or the spectrum of a given 
object, a comparison with our models might be useful as a de- 
scription of the overall ionization state and some bulk properties 
of the remnant, but it should be done with care. Moreover, the 
effect that the Raleigh-Taylor instabilities acting on the contact 
discontinuity (CD) might have on the ionization and heating of 
the outer ejecta layers is difficult to ascertain. As explained 
in Wang and ChevaUer (2001), the averaged density profile is 
smoothed by the instabilities, but the fingers of shocked ejecta 
retain a higher density than the surrounding shocked AM, so 
the results of our one dimensional simulations discussed in the 
following sections might not deviate much from a more realistic 
case. 

Even if spherical symmetry is assumed, it is worthwhile to 
note that, while SNIa progenitor systems are not supposed to 
substantially modify the surrounding ISM, the constant density 
hypothesis does not agree with current presupemova evolution 
models. These models imply a substantial modification of the 
AM, with profound effects in the dynamic evolution of the rem- 
nant, as discussed in Badenes and Bravo (2001). This possibil- 
ity has not been taken into account here, but will be explored in 
future work. 



Another issue that is of special concern is the validity of 
the adiabatic hypothesis in the hydrodynamic calculations: 
whereas radiative losses are generally not important for so- 
lar composition plasmas within time scales of a few thousand 
years, heavy element plasmas will radiate at a much faster rate, 
and the losses might have noticeable effects on the dynam- 
ics much earlier in the evolution of the SNR. The problem is 
complicated by the fact that losses are composition and ion- 
ization state dependent, so they can only be evaluated a pos- 
teriori, after the ionization calculations described in section 
3.3 are completed, and it is not possible to include their ef- 
fect in the hydrodynamics with our simulation scheme due to 
the fact that both calculations are decoupled. This a posteri- 
ori monitoring of the radiative losses in a heavy element NEI 
plasma can be done using the atomic data from Summers and 
McWhirter (1978) as described by Laming (2001). For eval- 
uation purposes, we define the time scale for the onset of ra- 
diative losses trad , as the time when the calculated a posteri- 
ori losses exceed 10% of the internal energy in a number of 
layers that amount to 5% of the total ejecta mass. This time 
scale is below 5000 years only for the deflagration models: 
trad^DEFa = 3.03 X 10^° s; t.ad.DEFf = 2.37 X 10^° s foT Wli inter- 
action with pam = 10"^"* g cm"-', with t^ad increasing for lower 
AM densities. Because radiative coohng is a runaway process, 
the validity of these models beyond ttad is difficult to determine, 
but it is worthwhile to notice that the losses are usually confined 
to a very small volume and the effect on the overall dynamics 
of the SNR should not be important (Hamilton, Sarazin, and 
Szymkowiak 1986). 

4. RESULTS 
4.1. Hydrodynamics 

The interaction of Type la ejecta with a constant density AM, 
based on a grid of one dimensional thermonuclear supernova 
explosion models, was first explored in Dwarkadas and Cheva- 
lier (1998, hereafter DC98). In that paper, the detailed dynam- 
ics of the interaction of six ejecta density profiles from ther- 
monuclear supernova models was examined and compared with 
three analytical density distributions: an exponential, a power 
law of index n = 7 with a constant density core, and a constant 
density profile. Except for the exponential, these analytical 
functions have been widely used before to approximate SNIa 
ejecta density as a function of velocity. Approximate temper- 
ature profiles were also calculated assuming solar abundances 
for the supernova ejecta. DC98 conclude that in all cases the 
density of the shocked ejecta increases from the reverse shock 
towards the contact discontinuity, that this density rise is usu- 
aUy coupled with a drop in temperature in the area close to the 
CD and that the sharp structures in the ejecta profiles, especially 
in the He detonation models, give rise to secondary waves prop- 
agating in the interaction region which could affect the instan- 
taneous X-ray emissivity of the remnants. 

When analyzing the ejecta-AM interactions, it is important 
to note that, once the ejected mass Me^, kinetic energy E]^ and 
density profile of an explosion model are fixed, the interaction 
with the AM follows a scahng law for the AM density pam 
(Gull 1973). We have tested this scaling law by computing the 
interaction of the supernova models with the three different AM 
densities mentioned in section 3.2. Our results follow the scal- 
ing law up to the precision of the hydrodynamic calculations. 
Therefore, it is sufficient to present the hydrodynamic calcula- 
tions for each model with a certain value of pam and then use 
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the characteristic magnitudes defined by DC98 (eqs. [3], [4], 
and [5]) to scale the results to any other pam value that might 
be of interest. 

^'=(#^1 . (3) 

t'=- = ^ (5) 

In Figure A3 we show our results for the ejecta-AM interac- 
tion with Pam = 10"^'*gcm"^. We have plotted the time evo- 
lution of the forward and reverse shock radii {rtwd, >'tev), the 
velocity of the forward shock Ufwd, the velocity of the reverse 
shock in the rest frame of the expanding ejecta Urev = (''rev/0 - 
{dr^ev/dt) and the expansion parameters for both shocks, de- 
fined as r]iwd,rev = d\n{riwa,rev) / d\n{t). The time axis spans be- 
tween 50 and 5000 years after the explosion, and the reverse 
shock parameters have been plotted only up to the time of re- 
bound. In fact, our simulations of the interaction of the ejecta 
with the AM start at 10^ s after the explosion, but we present 
our results only from year 50 on because,generaUy speaking, 
younger remnants are not expected to emit appreciably in X- 
rays. 

The dynamics of the forward shock is affected by the dif- 
ferences in the density profiles during the first thousand years, 
then all the models converge towards the Sedov-Taylor solu- 
tion (^7f^^.J = 0.4); this transition leads to the change of slope of 
the forward shock radii in the log-log plot (Fig. A3a). The 
shock trajectories of the deflagration models (DEFa, DEFf) lag 
behind the others because their is lower, and they can also 
be easily distinguished by their high r/ft^d values at early times, 
(Fig. A3e), about 50% higher than in the other models. Sud- 
den increases in rjf„d can be seen in the PDD models around 
/ = 3 X 10^ s, and in the SCH model at / = 6 x 10** s (the rise is 
outside the range of Figure A3e, but its effect can be seen at the 
beginning of the plot) and f = 6 x 10^ s. In general, high rji„d 
values are found in models that have high density material in 
the outermost ejecta, as the DBF models, or high density lay- 
ers in the outer ejecta preceded by a low density tail, as PDDa 
at r = 1.6 X 10^^ cm, PDDe at r = 1.4 x 10^^ cm, and SCH at 
r = 1.9 X let'em and r = 1.3 x 10^'' cm (see Fig. A2). These 
high density layers in the outer ejecta transfer their momentum 
to the shocked material and to the forward shock, leading to the 
increased r^f^d that we have noted. The high density material, 
once shocked, stays close to the contact discontinuity, and the 
low density tail, if present, is also recompressed by the reflected 
shocks that ensue when the high density layers are overcome 
by the reverse shock. The result is that the density enhance- 
ment effect close to the contact discontinuity is stronger for the 
DEF, PDD and SCH models than for the DDT and DET models, 
where the momentum transfer is more gradual (for examples of 
density maps see Fig. 3 in Badenes and Bravo 2003b). The 
dynamics of some of these models are compared to those of the 
power law and exponential analytical profiles in Badenes and 
Bravo (2003a). 

Summarizing, our conclusions agree almost completely with 
those of DC98: the density of the ejecta always peaks towards 
the contact discontinuity for all models at all times, and the rich 



internal structure of the explosion models produces a series of 
secondary shock waves that travel along the shocked ejecta and 
AM reheating and recompressing the material. Our results also 
agree with DC98 in that the mean temperature profile drops 
towards the contact discontinuity, but we have found that this 
behavior of the mean temperature cannot be extrapolated to the 
electron temperature, as is explained in the following section 
(see also Fig. 2 in Badenes and Bravo 2003b). We have also 
found that the models with dense layers in the outer ejecta, that 
is, the deflagration models and, to a lesser extent, the PDD mod- 
els and the SCH model, tend to achieve higher densities in the 
region behind the contact discontinuity. This wiU have a pro- 
found impact in the ionization and temperature calculations. 

4.2. Temperature, ionization, and emission measure 

The spectral characterization of a young SNR is a complex 
issue: since the density, ionization state, temperature and com- 
position of each fluid element are different, each region of the 
SNR will contribute differently to the total integrated spectrum, 
and so will each chemical element. A convenient way to mea- 
sure these contributions is the emission measure (EM) for ele- 
ment X, defined as 

EMx= I nxn^dV, (6) 

where Vsh is the volume of shocked ejecta. For identical physi- 
cal conditions and a common history, elements with equal emis- 
sion measures contribute equally to the total integrated ejecta 
spectrum. But the electron temperature and ionization time 
scale of the plasma, t = J n^dt, which play a key role in ther- 
mal NET spectra, are different for each fluid element, resulting 
in different spectra produced by fluid elements with identical 
emission measures. Borkowski, Lyerly, and Reynolds (2001) 
approached this problem for the shocked AM in Sedov SNRs 
by introducing distribution functions, plots of 7^ and r versus 
EM, but the use of this approach for the shocked ejecta would 
call for an individual distribution function for each chemical 
element due to the nonuniform chemical composition. An in- 
complete, yet meaningful, description can be achieved by tak- 
ing the first moment of the distribution functions and calcu- 
late, for each element X, an emission measure averaged elec- 
tron temperature (Tg)^ and ionization timescale (t)^. (Another 
quantity, an ionization timescale averaged electron temperature 
is generally necessary for a reliable modeling of X-ray spec- 
tra, but it is less important than {Ts)x and (t)^ ). We stress 
that we do not recommend single temperature, single ionization 
timescale spectral models to fit shocked ejecta in young SNRs, 
we are merely using the averaged quantities to describe average 
physical conditions of various layers of ejecta in our models. In 
all the spectra and calculations presented in the following sec- 
tions we have used the full profile of electron temperatures and 
ionization timescales throughout the ejecta, not the averaged 
quantities. 

The EMx evolution is shown in Figure A4 for pam = 
10~^"*gcm"^, between 20 and 5000 years after the explosion. 
The contributions from different chemical elements to the 
ejecta spectra depend strongly on both the composition profile 
of the models and their dynamic evolution. The density en- 
hancement effect towards the contact discontinuity makes the 
chemical elements in the outer layers of the ejecta more promi- 
nent than those in the inner layers, so the importance of Fe in 
the model spectra is generally much less than might be expected 
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on the basis of a Type la elemental composition alone. This 
relatively low prominence of Fe in spectra of Type la SNR can- 
didates might have been noticed on a number of occasions, of- 
ten accompanied by inordinately high apparent abundances of 
other elements (see Hendrick, Borkowski, and Reynolds 2003 
for the SNRs 0548-70.4 and 0534-69.9 and Lewis et al. 2003 
and van der Heyden et al. 2002 for N103B). The estimated 
ejected masses of various elements, which often indirectly rely 
on the assumption that chemical abundances are proportional to 
the fitted EM for each element in the spectrum, are difficult to 
reconcile with the yields of theoretical explosion simulations. 
This discrepancy between EM and ejected mass is most dra- 
matic in the deflagration models (Fig. A4 e and f), whose spec- 
trum is completely dominated by C and O, with emission mea- 
sures of Fe about two orders of magnitude lower at all times, 
even though the ejected mass of Fe is higher than that of C or 
O. Note, however, that the peak value of the EM-pe in the defla- 
gration models is about the same as in the other models, with 
the exception of PDDa, and that radiative cooling could reduce 
considerably the EM of C and O in the deflagration models 
(see section 3.5). The prompt detonation model (Fig. A4g) is 
the only one whose spectrum is clearly dominated by Fe at all 
times, while in DDTa and PDDa (Fig. A4a and c) Fe takes over 
only after a few hundred years, with important contributions of 
Si and S (and O for DDTa) throughout the SNR evolution. The 
rise in the emission measures around ? = lO'^ s is due to the 
propagation of the reverse shock after it rebounds at the cen- 
ter, reheating and recompressing the ejecta; this rise happens 
earlier and is more gradual for the elements in the inner layers 
than for those in the outer layers. Model SCH is dominated 
by Fe only when the reverse shock is propagating through the 
He detonation layer, afterwards O takes over, because a second 
density enhancement region forms behind the interface between 
this layer and the rest of the exploded WD. 

Since the ionization and electron heating processes proceed 
faster at higher densities, the corresponding (7i)x ^^"1 (r)^ 
plots in figures A5 and A6 are also affected by the enhance- 
ment towards the contact discontinuity and the reverse shock 
rebound, albeit in different ways. Those models with stronger 
density enhancement close to the CD (DEFa, DEFf, PDDe, and, 
to a lesser extent, PDDa and SCH) tend to have higher ion- 
ization timescales for most elements, and in general, elements 
closer to the CD have higher (r)^. In the deflagration mod- 
els, for instance, C and O are always at a higher r and initially 
hotter than other elements. The {T^)q q drops at later times be- 
cause in these models the density is so high that electron-ion 
temperature equiUbration is achieved for most of the C and O 
in the ejecta before t = lO'" s and afterwards the electrons just 
cool due to adiabatic expansion of this region of the SNR. Iron 
is generally hotter and at a higher t in the DDT and PDD mod- 
els than in the DEF models. The anomalous behavior of the 
plots for some elements (for instance, Ni in DDTa and PDDa 
or Ca in SCH) is due to the averaging in EM and can be un- 
derstood by comparing the curves with the chemical compo- 
sition profiles of Fig.Al. As the reverse shock advances into 
regions with a much higher concentration of a given element, 
the newly shocked (and therefore cooler and less ionized) lay- 
ers soon dominate the EMx and the averaged quantities shift 
their values accordingly. It is worth noting that the enhanced 
electron heating rate due to the higher densities towards the CD 
compensates for the lower specific internal energies found in 
that region, so we observe 7^ profiles that always peak at the 



CD (i.e. electrons and ions are closer to thermal equilibrium at 
the CD, whereas <C 7] behind the reverse shock) in contrast 
to the behavior of the mean temperatures (DC98). 

Our results for lower AM densities are similar, and we have 
found that the hydrodynamical scaling laws (section 4.1) also 
work reasonably well for the emission measures, except at ex- 
tremely low AM densities. Scaled EMx agree to within a factor 
of 2, and so do the {t)^. The temperature scaling seems more 
complex. The deviations are accounted for by the difference in 
average ionization stale which is expected in SNRs that evolve 
in different pam, because hydrodynamical scahng does not ap- 
ply to ionization and electron heating processes. A detailed 
analysis of the ambient density effect on the spectral character- 
istics of our models will be the subject of future work. 

4.3. Synthetic spectra 

The integrated synthetic spectra from the ejecta in our SNR 
models are presented in figure A7 for pam = 10~^* g cm"^, 500, 
1000, 2000, and 5000 years after the explosion (without inter- 
stellar absorption). As explained in the previous section, the 
contributions from the different elements to the spectra depend 
strongly on the details of the hydrodynamic evolution, and in 
particular on the density enhancement towards the CD. Note 
how the spectra of the deflagration models are always domi- 
nated by C and O, to the point that their continua 'veil' the 
lines of the other elements at early times. Even though their 
composition profiles are similar, PDDa and PDDe have a richer 
line spectrum than DDTa and DDTe respectively, because their 
mean ionization state is more advanced due to the higher den- 
sities close to the CD (section 4.1). The strength of the Fe Ka 
line varies from model to model, being important at all times in 
DDTa, PDDa, DET and SCH, which have more Fe in flie outer 
layers of ejecta, while for DDTe and PDDe Fe Ka is noticeable 
only after the first 1000 years, for DEFf after 2000 years and 
much later for DEFa. Oxygen can be easily identified in the 
models with high EMq and low EMpe (otherwise its presence 
is partiaUy masked by the Fe L complex): DDTe and PDDe at 
early times, SCH up to 2000 years and the deflagration models 
at all times. These oxygen bright models with weak Fe lines 
could be easily mistaken for core-collapse SNRs. For the de- 
flagration models in particular, their strong C and O continua 
could make spectral analysis difficult in presence of substantial 
interstellar absorption and in view of calibration problems and 
low spectral resolution of CCD detectors at low photon ener- 
gies. 

4.4. Comparison with observations: the Tycho SNR 

A comparison with observations is essential for assessing 
validity of theoretical models as a spectral characterization of 
young type la SNRs and for learning about their progenitors. 
Our models are well suited for qualitative comparison with ob- 
served spectra, but a quantitative, detailed study of a given ob- 
ject is much more difficult. The exceUent quality of present 
day observations from XMM-Newton and Chandra has sur- 
passed our ability to model them accurately even with sophis- 
ticated theoretical models, and additional issues such as mod- 
eUng of shocked AM spectrum, nonthermal spectrum subtrac- 
tion, matching shock velocities, expansion parameters and ap- 
parent radii, and spatially resolved spectroscopy of the ejecta 
will have to be considered. Analysis of each particular SNR 
might call for some modifications in the calculations or for a 
more exhaustive exploration of the parameter space, and this 
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task is deferred to future work. However, a preliminary com- 
parison of our results with observations of a well studied SNR 
is very useful, as it allows us to assess the strengths and Umita- 
tions of our models. 

The best object for this comparison is Tycho, the remnant 
of SN 1572. Tycho is considered the remnant of a Type la 
supernova, and it has been extensively studied in radio, opti- 
cal, and X-rays. The X-ray emission from the ejecta has been 
observed with high angular and spectral resolution by XMM- 
Newton and Chandra (Decourchelle et al. 2001; Hwang et al. 
2002), and spatially integrated but with moderate spectral res- 
olution by ASCA (Hwang, Hughes, and Petre 1998; Hwang 
and Gotthelf 1997, hereafter HG97). However, these extensive 
observational studies of Tycho resulted in just a few compar- 
isons with SNIa explosion models. Itoh, Masai, and Nomoto 
(1988), and Brinkmann et al. (1989) calculated the predicted 
X-ray spectra at the age of Tycho using the deflagration model 
W7 (Nomoto, Thielemann, and Yokoi 1984), and compared it 
to the Tenma and EXOSAT data to find that the observations 
could not be reproduced by the W7 model unless substantial 
modifications to its structure where introduced. 

We have made no attempt to reproduce the apparent size of 
the remnant and the observed proper motion of Tycho; we just 
compared the X-ray spectra of some of our models with the spa- 
tially integrated XMM spectrum (Decourchelle et al. 2001, Fig- 
ure 8). The contribution from the shocked AM has been mod- 
eled with the Sedov model in XSPEC (Borkowski, Lyerly, and 
Reynolds 2001), with mean shock temperatures and shocked 
ISM emission measures derived from our hydrodynamical cal- 
culations for each case (see Tab. A2), and assuming no colU- 
sionless electron heating at the shock (kTe; = 0), in agreement 
with optical spectroscopy (Ghavamian et al. 2001). The ioniza- 
tion age, defined as the product of postshock electron density 
and the remnant's age, is 2.93 x 10^'^scm"^. We note that the 
shock speeds we obtain (and therefore the mean shock tempera- 
tures kT^h used in the Sedov models) are compatible with X-ray 
measurements (Hughes et al. 2000), but not with optical or ra- 
dio observations (Hwang et al. 2002, and references therein). 
The distance to Tycho was assumed to be 2.3 kpc, and we have 
used a A^h of 0.45 x lO^^cm'^ as in HG97. 

A simple visual comparison of Tycho's spectrum with the 
model spectra at t=500 yr (Fig. A7) reveals that one can ehm- 
inate the DET model (absence of Si, S, and Ca lines) as well 
as DEFa, DEFf and DDTe (mainly due to the absence of the 
FeKa line). The contribution from the ISM (Sedov models) to 
the ejecta spectra varies a lot depending on the type of explo- 
sion. It is quite large, specially for the high energy continuum, 
in models DDTa, DDTe and DET, and quite small in the others, 
except PDDa which is an intermediate case. For the DEF mod- 
els in particular, almost all the emission comes from the ejecta, 
while for DET and DDTe the contribution to the underpredicted 
Unes from the ISM is small, and therefore these models remain 
unable to account for the spectrum of Tycho. The spectra of 
the most promising models have been plotted alongside the ob- 
served spectrum in Figure A8. 

For a more quantitative comparison we have computed the 
main characteristics of the more prominent emission Unes of 
our models, and compare them to the observational results of 
HG97 in Tables A2 and A3. The flux for the Ka complexes of 
Si, S, Ca and Fe was obtained by adding the contributions from 
all lines due to atomic transitions from n = 2 to n = 1 . 

We concentrate now on models DDTa, PDDa, PDDe, and 



SCH, which have not been eliminated previously. PDDa, 
PDDe, and SCH give too strong Si and S Lya Unes, with re- 
spect to the corresponding Ka Une complexes. TUese three 
models have in common the formation of an extended high- 
velocity envelope (corresponding to the accreted He envelope 
in the SCH model and to the external layers detached from the 
white dwarf in the pulse previous to the formation of the det- 
onation in the PDDa and PDDe models), which is separated 
from the rest of the ejecta by an abrupt change in density and 
followed by dense layers rich in Si and S (see Fig. Al). These 
dense layers remain dense when shocked, speeding up the ion- 
ization processes inside them and producing more highly ion- 
ized Si and S and therefore higher Lya/Ka ratios. Other model 
ratios compare reasonably well with the observed values, most 
of them within a factor two to three. The Fe Ka line centroids 
sUow deviations larger than 50 eV only for models PDDa and 
SCH. 

Given the limitations of the present work and the consider- 
able uncertainties in the atomic data it is impossible to reach 
a definitive conclusion, but our preliminary results point to a 
delayed detonation explosion as the progenitor of Tycho, prob- 
ably similar to the DDTa model, with a smooth density profile 
and a large amount of Fe in the outer ejecta. Pure deflagration 
and pure detonation models can be confidently eUminated. The 
only features that model DDTa does not reproduce weU are the 
Ca Ka / Si Ka ratio and the S K/3 / S Ka ratio, suggesting that 
the ionization state of S in Tycho may be a little bit higher than 
in the model. The angular size of Tycho (8') implies a radius 
of 2.8(Z)/2.3 kpc) pc, while the radius of the DDTa model at the 
age of Tycho is about 3.2 pc (see Fig. A3). If the distance 
has been estimated correctly, then p^m must be higher than 
10~^'*gcm""'. In simulations with a moderately higher value 
of pam, the ionization state of S in the DDTa model increases, 
bringing the S K/3 / S Ka ratio closer to the observed value 
while keeping the rest of the S and Si line ratios within reason- 
able limits. A more detailed study of Tycho including these and 
other issues will be the subject of future work. 

5. CONCLUSIONS 

We have performed self-consistent hydrodynamic and ion- 
ization state simulations of the interaction of a grid of SNIa ex- 
plosion models with a uniform AM, producing a set of synthetic 
spectra for the integrated thermal X-ray ennission from the 
shocked ejecta in young SNRs. The spectra show remarkable 
differences between the models: line strengths and strength ra- 
tios, overall luminosity, shape of the continuum and other pa- 
rameters depend strongly on the type of supernova explosion 
that gave birth to the SNR. These differences stem from the hy- 
drodynamic evolution, that is determined by the density profile 
of the ejecta synthesized in the explosion, and from the evo- 
lution of the electron heating and plasma ionization processes, 
which depend on both the dynamic history and the chemical 
composition of the ejecta. A close cormection is thus estab- 
lished between supernova explosions and young SNRs, im- 
plying that SNRs have the potential to become a new tool to 
discriminate among various SNIa explosion models, a possi- 
biUty that should be tested by comparing models with obser- 
vations. It is also clear that any analysis of a thermal X-ray 
spectrum from the shocked ejecta in a young SNR, whatever 
its type, has to make reaUstic assumptions about the density 
and chemical composition profile for the ejecta, because us- 
ing uniform densities or homogeneous compositions is unreal- 
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istic and will produce misleading results. The density enhance- 
ment effects caused by the hydrodynamic evolution can lead 
to emission measure ratios very different from elemental ratios 
(by number) in the material ejected by the explosion. This is 
because the density enhancement amplifies the contribution to 
the emitted spectrum by the elements in the outermost layers 
of ejecta, where the differences between the models are most 
pronounced. In contrast with the results of Dwarkadas and 
Chevalier (1998) for the mean gas temperature, we find that 
both density and electron temperature are highest at the contact 
discontinuity in all cases. 

We have made a preliminary comparison between our theo- 
retical spectra and the observed spectrum of the Tycho SNR. 
The best choice for Tycho is a delayed detonation model with 
a large kinetic energy and a high iron content. In particu- 
lar, our delayed detonation model DDTa provides reasonable 
agreement with observed line strengths, Une ratios, and Fe Ka 
line centroid. Explosion models characterized by high den- 
sity layers in the outer ejecta preceded by an extended low- 
density envelope, like pulsating delayed detonations and sub- 
Chandrasekhar models give rise to an excess of ionization of 
Si and S, which leads to high Lya/Ka line ratios. Deflagration 
models characterized by a large buffer of C and O in the exter- 
nal layers would produce spectra with too weak lines and too 
strong continua. One important result is that the emission mea- 
sure ratios (as well as the line ratios) cannot be safely used as 
a measure of the relative abundances of various elements in the 
supernova ejecta; they can in fact differ by as much as two or 
three orders of magnitude from the actual abundance ratios. 

A potentially important issue is related to various hydrody- 
namical instabilities occurring during the explosion itself or 
shortly thereafter. Multi-dimensional simulations of Type la 
explosions show the presence of these instabilities (Khokhlov 



2000; Gamezo et al. 2003; Niemeyer et al. 2003; Garcia-Senz 
and Bravo 2003; Bravo and Garcia-Senz 2003), which may lead 
to inhomogeneities in freely expanding ejecta. There is indeed 
observational evidence for the presence of fast-moving ejecta 
clumps in Tycho SNR (Hwang and Gotthelf 1997). But a quan- 
titative understanding of the dumpiness of the freely expanding 
ejecta is lacking at this time, both from theoretical and obser- 
vational perspective (see, for instance, Gamezo et al. 2003, and 
Thomas et al. 2002). A small degree of dumpiness is not likely 
to affect our results based on 1-D simulations, but the situation 
will be different for highly inhomogeneous ejecta. In extreme 
cases, dense clumps of shocked ejecta might significantly con- 
tribute to spatially-integrated X-ray spectra. The problem of 
clumpy supernova ejecta is clearly outside of the scope of this 
work, and should be addressed in the future through multidi- 
mensional simulations and through detailed observational stud- 
ies of Type la SNR morphologies. 

All the models, plots, and spectra shown here are available 
from the authors upon request. 
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APPENDIX 

THE SUPERNOVA EXPLOSION MODELS 

In this appendix we give more details about the SNla explosion models. As has been mentioned before, the hydrocode, the 
nucleosynthesis, and the light curve codes, as well as the physics included (equation of state, nuclear reaction rates, etc) are the same 
as described by Bravo et al. (1996). The parameters of the models are given in Table Al. 

The detonation model was obtained from a ~ 1.38M0 WD in hydrostatic equilibrium (composed of equal masses of ^■^C and ^^O 
plus a 1 % by mass of ^^Ne) whose internal energy structure was adjusted to an adiabatic thermal gradient. The ignition was initiated 
by incinerating the mass in the central layer, and afterwards the detonation propagation was obtained consistently by solving the 
hydrodynamic and nuclear evolutionary equations. Details of a similar model can be found in Bravo et al. (1996). 

The sub-Chandrasekhar mass model was obtained from a WD formed by a C-O core of 0.8 M© surrounded by a He envelope 
of O.2M0. This envelope was the result of He accretion over the C-O core at a steady rate of 3.5 x lO'^M© yr"'. The hydrostatic 
evolution of the white dwarf subject to accretion was computed by J. Jose, who kindly provided us with the initial model for the 
supernova explosion calculation (private communication). In this initial model, the temperature at the base of the He envelope was 
high enough to induce a spontaneous He detonation. The evolution past this point was followed with the same hydrocode as above. 

The rest of the explosion models started from the same initial configuration: an isothermal white dwarf in hydrostatic equilibrium, 
with the same chemical composition and central density as that used in DET model (the differences in the thermal structures of both 
configurations account for the small difference in total mass that can be seen in Table Al). Explosion models starting from different 
central densities pc (i.e. WD masses) do not produce substantially different energies or Ught curves (with the exception of a slight 
decrease of ^^Ni yield with increasing central density due to a larger electron capture rate. Bravo et. al. 1993), so we did not consider 
variations in p^. 

The flame propagation velocity in the pure deflagration models was obtained as the maximum between the laminar flame velocity 
(as given by Timmes and Woosley 1992 and updated by Bravo and Garcia-Senz 1999) and the turbulent velocity, vrj. The turbulent 
velocity was calculated as vr?- = nru/r^T, where rg is the flame radius, trj is the local Rayleigh-Taylor time scale at the flame 
location, and k is a parameter given in Table Al (see Bravo et. al. 1996 for details). Furthermore, the mixing of matter and energy 
across the flame front was limited as proposed by Wheeler et al. (1987). In this formulation, the mixing is not allowed until the 
Rayleigh-Taylor front has propagated to a prescribed fraction, 9, of the mass of the shell (in all our deflagration models we took 9 = 
0.5). Thus, if the mass of the sheU ahead of the flame front is AM, the mixing was started when the condition ^^^dt = 6x AM was 
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held, with to equal to the time of incineration of the last incinerated shell, and with = 47rr|, x ptv x vr^. Once this condition was 
fulfilled, the transfer of internal energy between the Rayleigh-Taylor unstable shells was allowed. Afterwards, the flame propagation 
was obtained consistently by the consequent increase in the nuclear energy generation rate and its feedback on temperature. At 
densities below a few times 10^ g cm~^ the energy generated goes predominantly to create and maintain an electron-positron pair gas 
rather than to increase the temperature, which is the ultimate cause of the quenching of the flame. 

In the delayed detonation models the flame was propagated initially as a deflagration. In that phase, the flame velocity was taken 
as a constant fraction, t, of the local sound velocity (in our models we took i = 0.03). As the flame traveled through lower and lower 
densities, the sound velocity and hence the deflagration speed decreased, and the flame was eventually quenched. Afterwards, the 
flame front was artificially accelerated up to a large fraction of the sound speed. Following this fast propagation through a few mass 
shells (typically 2-3 shells were enough) a detonation formed and propagated through the rest of the star. The location of the transition 
is univocally determined by the parameter ptr (Tab. Al), that is the density ahead of the flame at which the sudden acceleration was 
imposed. In delayed detonation models as well as in the pure detonation and sub-Chandrasekhar models, Rayleigh-Taylor mixing 
was not allowed. The algorithm for the flame propagation in pulsating delayed detonation models was the same as in delayed 
detonation models, the only difference being that in PDD the detonation was not triggered until the WD had pulsated. Further details 
(in particular the nucleosynthetic output and more complete information about the light curves) are available to interested readers 
upon request. 
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Table A 1 

Properties of the explosion models 



Model 


Parameter^ 










Mc+o 


Ms, 


Ms 


Ma, 




Mf, 






(Mo) 


(10" erg) 


(mag) 


(mag) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


DET 




1.38 


1.59 


-19.87 


1.01 


0.0003 


0.0008 


0.0007 


0.0002 


0.0002 


1.22 


SCH 




0.97 


1.01 


-17.53 


1.03 


0.09 


0.14 


0.10 


0.024 


0.026 


0.50 


DEFa 


0.06 


1.37 


0.51 


-18.97 


0.85 


0.69 


0.025 


0.017 


0.0038 


0.0040 


0.55 


DEFf 


0.16 


1.37 


0.84 


-19.43 


1.04 


0.48 


0.017 


0.012 


0.0028 


0.0029 


0.75 


DDTa 


3.9 X 10' 


1.37 


1.40 


-19.73 


1.11 


0.04 


0.087 


0.071 


0.019 


0.022 


1.03 


DDTe 


1.3 X 10' 


1.37 


1.02 


-19.00 


0.94 


0.19 


0.25 


0.19 


0.046 


0.054 


0.56 


PDDa 


4.4 X 10' 


1.37 


1.45 


-19.79 


1.10 


0.02 


0.055 


0.045 


0.012 


0.015 


1.11 


PDDe 


7.7 X 10* 


1.37 


1.12 


-19.02 


0.95 


0.10 


0.27 


0.22 


0.057 


0.067 


0.58 



'The parameter given is k for DEF models and pt, for DDT and PDD models (see the appendix for definitions) 
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Table A2 

MODhL SPHC I RAL hHAl LiRhS AT I IIH AGH Oh TYCIIO SNR 



Model 


Si Ka line flux 


Fe Ka line centroid' 




(£Mis„/47rD2' 




(lO^'phcm-^s'') 


(keV) 


(keV) 


(cm-5) 




52.7 


6.458 






DDTa 


11.95 


6.436 


27.4 


1.39 X 10" 


DDTe 


22.60 




24.0 


1.16 X 10'' 


PDDa 


23.22 


6.537 


28.4 


1.45 X 10" 


PDDe 


90.50 


6.465 


26.8 


1.22 X 10" 


DEFa 


3.46 




22.6 


0.58 X 10" 


DEFf 


21.02 




25.7 


0.93 X 10" 


DET 


1.89 


6.438 


27.9 


1.48 X 10" 


SCH 


31.35 


6.510 


25.1 


1.21 X lO" 



^Centroids marked as — would be impossible to determine from the model spectra due to 
the weakness of the line 

''Parameters for the Sedov model for the shocked ISM 
•^Observational results of Hwang and Gotthelf (1997) 



Table A3 

Diagnostic line ratios for the models compared with the results of HG 97 



Line Ratio 


SiK(/3 + 7)/ 


Si Lya/ 


SK/3/ 


SLya/ 


SKa/ 


Ca Ka/ 


FeKa/ 




SiKa 


SiKa 


SKa 


SKa 


SiKa 


SiKa 


SiKa 


HG97 


0.129 


0.028 


0.065 


< 0.010 


0.26 


0.010 


0.008 


DDTa 


0.090 


0.013 


0.011 


0.005 


0.55 


0.045 


0.006 


DDTe 


0.105 


0.012 


0.012 


0.003 


0.33 


0.014 


0.001 


PDDa 


0.161 


0.143 


0.068 


0.098 


0.49 


0.086 


0.029 


PDDe 


0.162 


0.215 


0.088 


0.079 


0.50 


0.042 


0.001 


DEFa 










0.34 






DEFf 


0.147 


0.146 






0.27 


0.028 




DET 


0.069 


0.026 






0.49 




0.074 


SCH 


0.133 


0.175 


0.070 


0.050 


0.29 


0.007 


0.005 



Note. — Ratios marked as — would be impossible to determine from the model spectra due to the 
weakness of the lines involved 



FIG. A 1 .— Chemical composition profile of the ejecta for the models DDTa (a), DDTe (b), PDDa (c), PDDe (d), DEFa (e), DEFf (f), DET (g), and 
SCH (h). Elements are represented as follows: C is the thin dash-triple-dotted plot (green in the electronic edition), O is thick dash-triple-dotted 
(blue), Si is thick short-dashed (cyan), S is thin short-dashed (yellow), Ar is long-dashed (purple), Ca is dash-single-dotted (gray), Fe is thick solid 
(red) and Ni is dotted (brown). The thin solid plot (black dotted plot in the electronic edition) is the normalized density, p/pmax- 



Fig. A2. — Density profiles for the supernova models lO' s after the explosion. Model DDTa is the thin solid line, DDTe is dotted, PDDa is thin 
short-dashed, PDDe is dash-single-dotted, DEFa is dash-triple-dotted, DEFf is long-dashed, DET is thick short-dashed, and SCH is thick solid. 

Fig. A3. — Temporal evolution of the shock parameters rt„d (a), (b), ut„d (c), Mrev (d), r)!„d (e) and r}rev (f) for an interaction with pam = 
lO"^'' gcm~^. Models marked as in Figure A2. 

FIG. A4.— EMx/4ttD^ evolution for the models DDTa (a), DDTe (b), PDDa (c), PDDe (d), DEFa (e), DEFf (f), DET (g), and SCH (h), for 
a distance of D=10kpc. The thin solid plot (black in the electronic edition) corresponds to the total EM/AitD^ of the ejecta (i.e., their total 
brightness). Elements labeled as in Figure Al. C and O plots overlap in the deflagration models. 

FIG. A5.— Emission measure averaged ((Ti)^) for the models DDTa (a), DDTe (b), PDDa (c), PDDe (d), DEFa (e), DEFf (f), DET (g), and 
SCH (h). Elements labeled is as in Figure Al. C and O overlap in the deflagration models. 



Fig. A6.— Emission measure averaged r ((r)^)for the models DDTa (a), DDTe (b), PDDa (c), PDDe (d), DEFa (e), DEFf (f), DET (g), and 
SCH (h). Elements labeled as in Figure Al. C and O overlap in the deflagration models. 

Fig. A7.— Spatially integrated synthetic ejecta spectra for the models DDTa (a), DDTe (b), PDDa (c), PDDe (d), DEFa (e), DEFf (f), DET (g), 
and SCH (h). Spectra plotted with a solid line correspond to 500 years after the explosion (blue in the electronic edition), dotted line to 1000 years 
(red), dashed line to 2000 years (orange) and dash-dotted line to 5000 years (green). Fluxes are for objects at a distance of £)=10kpc. Spectra are 
convolved with the response of the EPIC-MOSl camera onboard the XMM-Newton satellite. The Ka lines of Fe, Ca, S, and Si, as well as the O 
Lya line, have been marked for clarity. 



X-rays from thermonuclear SNRs 



11 



Fig. A8. — Top: Spectrum of Tycho observed by the XMM EPIC-MOSl camera (Decourchelle et al. 2001, reproduced with permission). Bottom: 
Predicted spectra for the models DDTa, PDDa, PDDe and SCH. The solid Une is the total spectrum, the dotted line is the ejecta spectrum, and the 
dashed hne is the ISM spectrum. 
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